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The ultrasonic wave propagation characteristics were measured for IN-100, a 
powder metallurgy alloy used for aircraft engine components. This material was 
selected as a model system for testing the feasibility of characterizing the micro- 
structure of a variety of inhomogeneous media including powder metals, ceramics, 
castings and . composites. The data were obtained for a frequency range from about 2 
to 20 MHz and were statistically averaged over numerous volume elements of the sam- 
ples. Micrographical examination provided size and number distributions for grain 
and pore structure. The results showed that the predominant source for the ultra- 
sonic attenuation and backscatter was a dense (~ 100/mm 3 ) distribution of small 
micropores (~ 10 pm radius). Two samples with different micropore densities were 
studied in detail to test the feasibility of calculating from observed microstruc- 
tural parameters the frequency dependence of the microstructural backscatter in the 
regime for which the wavelength is much larger than the size of the individual scat- 
tering centers. Excellent agreement was found between predicted and observed values 
so as to demonstrate the feasibility of solving the forward problem. The results 
suggest a way towards the nondestructive detection and characterization of anomalous 
distributions of raicropores when conventional ultrasonic imaging is difficult. The 
findings are potentially significant toward the application of the early detection 
of porosity during the materials fabrication process and after manufacturing of 
potential sites for stress induced void coalescence leading to crack initiation and 
subsequent failure. 


INTRODUCTION 


The acquisition and interpretation of ultrasonic data for the inference of 
microstructural information is not new and has been carried out by several re- 
searchers in the recent past [ref. 1-7]. The objective of this work is to apply 
these techniques to the types of powder metallurgy alloys used in the manufacture of 
aircraft engine disc components. The particular application motivating this work is 
the development of an approach to the nondestructive detection and characterization 
of anomalous microstructure when conventional ultrasonic imaging is difficult. Such 
an approach could be important tracking porosity during the materials process as 


*This work was sponsored in part by Rockwell International Research and Development 
Funds and in part by the Center for Advanced Nondestructive Evaluation, operated by 
the Ames Laboratory, USDOE, for the Defense Advanced Research Projects Agency and 
the Air Force Wright Aeronautical Laboratories/Materials Laboratory under Contract 
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well as during service as a finished part in the early locating of potential sites 
for crack initiation. 

The material chosen for these tests is IN-100, a Ni-base alloy prepared by 
Pratt & Whitney Aircraft Co. This material is used in the F-15 and F-16 fighter 
aircraft engine at a rate of 1000 lb of billet material per engine to manufacture 
three compressor discs, four turbine discs and two spacers. The following sections 
discuss the experimental approach, the analytical models and the comparison between 
theory and experiment. 


SAMPLE CHARACTERIZATION 


This section treats the micrographical and ultrasonic characterization of the 
samples, the determination of bulk microsructure parameters from the 2D optical 
micrographs and he ultrasonic measurement technique [ref. 8]. 


Micrographical Characterization 

Suitable samples were prepared from a number of different parent samples sup- 
plied by Pratt & Whitney Co. Selected sample surfaces were polished and examined 
under the optical microscope. Figure 1 shows an optical micrograph of the grain 
structure. Figure 2 gives histogram data displaying the grain size distribution for 
two samples. These data show a uniform grain size of about 2-4 p,m. 
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Figure 1 

Photomicrograph of IN-100 grains. 


Figure 2 

Grain size distribution in IN-100. 


One of the samples tested (sample C) is a section of a used F100 engine 2nd 
disc which was pulled out for spin pit tests. The section is shown in cross-section 
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in figure 3 together with lines drawn to indicate traverses with the microscope to 
study surface features# Figure 4 shows sample microphotos taken along some of these 
traverses indicating the presence of micropores. The microporosity is a by-product 
of the powder metallurgy process and has its origin in hollow powder particles with 
argon gas trapped during the atomization process. Although a sparse, uniform dis- 
tribution of micropores is not thought to be the limiting factor in low cycle fa- 
tigue strength, the micropores do contribute significantly to the ultrasonic noise 
background and must be taken into account. Table 1 presents a summary of the micro- 
pore populations found along the various traverses. Statistical analysis of the 
micrograph data indicate that the pore distribution is not random but show a ten- 
dency towards clustering. Figure 5 presents the information of Table 1 in graphical 
manner. Figures 6a and 6b give pore size distribution histograms for two of the 
samples tested. Sample A was found to have an approximate raicropore density of 
330/mm^ whereas sample C had 90/mm^. For sample C the majority of micropores fall 
in the range of between 5 and 10 pm in diameter whereas in sample C the diameters 
range from ~ 1 pm to 50 pm . 
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Figure 3 Section of F100 engine disc. 


In order to employ the histogram data obtained from the optical examinations to 
estimate volume micros tructural information, the expressions for the mean volume 
porosity p v and the expectation value for the sixth power of the pore radius E(r 6 ) 
were derived: 
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Figure 4 Microphotos of area elements (along tranverses of figure 3) 
showing micro porosity. 
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Figure 5 Graph of number of pores per unit length, number of pores per unit 

area and pore separation plotted vs distance from disc center to run 
(as indicated by traverse number, see Fig. 3). 
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Figure 6 Pore size distribution (a) sample C - IN-100, 
(b) sample A - IN-100. 
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Table 1. 

Summary of 

Observations in Spin 

Test 

Sample 


Traverse 

Number 

Traverse 

Traverse Length 
Symbol mm 

Sample 

Length 

mm 

Number 

of 

Frames 

Sample 

Area 

(mm 2 ) 

Number 

of 

Pores 

Number 

of 

Pores/ 

Linear 

mm 

Number 

of 

Pores 

2 

mm^ 

Pore 
Spacing 
(Av) mm 

1 

K 

40 

36.9 

81 

21.9 

29 

0.79 

1.3 

1.27 

2 

J 

40 

39.6 

87 

23.5 

28 

0.7 

1.2 

1.4 

3 

I 

40 

40.6 

89 

24 

54 

1.3 

2.3 

0.75 

4 

H 

16 

15 

33 

8.9 

27 

1.8 

3.0 

0.55 

5 

A 

8 

7.3 

16 

4.3 

10 

1.4 

2.3 

0.7 

6 

B 

8 

8.2 

18 

4.9 

3 

0.36 

0.6 

2.7 

7 

C 

19 

18.2 

40 

10.8 

35 

1.9 

3.2 

0.5 

8 

D 

10 

10.5 

23 

6.2 

15 

1.4 

2.4 

0.7 

9 

E 

10 

10.5 

23 

6.2 

14 

1.3 

2.3 

0.75 

10 

F 

21 

16 

31 

9.5 

34 

2.1 

3.6 

0.47 

11 

G 

21 

15.5 

_34 

9.2 

22 _ 

1.4 

2.4 

0.7 

Total Sample 


218 


129 

271 

1.27 

0.87 
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where n2(r')dr' is the number of circles per unit area (with radii between r' and r' 
+ dr' as measured from micrographs. Each circle represents a cut through a pore by 
the polished sample surface. 113 (r*) is the number of micropores per unit volume 

element, Ar. is the bin width of the pore radius histogram, r. is the average pore 

radius for each bin and n^ is the number of pores in each bin. 

These formulae were used for calculations in the detailed examination of two 

demonstrative samples discussed below. 


Ultrasonic Characterization 

The samples were studied in two ways. The ultrasonic attenuation was measured 
in the frequency range from about 5 to 100 MHz on samples A and B with nominally 
similar pore size distributions. The shear wave backscattered power was measured in 
the frequency range 2—20 MHz on samples C and A described in detail in the text. 

In the ultrasonic attenuation measurements the data were obtained in the 10 MHz 
to 100 MHz range of frequencies in several smaller sub-bands. A low attenuation 
buffer (A1 or fused silica) was typically used to obtain the transducer response and 
the reflection coefficient of the sample interface. The technique and calculations 
are illustrated in figure 7 for the case of longitudinal waves. For the case of 
shear waves, an additional technique was employed for rapid data acquisition and is 
illustrated in figure 8. Here the buffer is water and the angle of incidence is 
such that the longitudinal waves in the water produce shear waves in the sample. 
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TRANSDUCER 
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Figure 7 Sample configuration for measurement of longitudinal wave attenuation. 
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Figure 8 Sample configuration for measurement of shear wave attenuation. 

Because of the possible generation of leaky waves generated at the interfaces and 
the non-circular beam spreading associated with non-normal incidence, the measure- 
ments were augmented with calibration measurements, carried out with narrow band 
quartz shear wave transducers in direct tightly bonded contact with the sample. 
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db/cm 


In reference to figure 8, Toi, Tiq are the water-sample, sample-water trans- 
mission coefficients; a g is the sample shear wave attenuation in nepers/cm, is 
the water longitudinal attenuation in nepers/cm, l g ^, l g 2 are t ^ ie s ^ ear wave paths 
in the sample under test, 1 l w 2 are the wave path length in water* 

The measurements are summarized in figures 9 and 10 which display attenuation 
values as a function of frequency for two different samples and for three different 
locations for the same sample. 
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Figure 9 Figure 10 

Longitudinal (L) and shear (S) wave Longitudinal wave attenuation at three 

attenuation for 2 samples of IN-100. different locations on sample C. 


Table 2 summarizes the information collected by given representative values for 
attenuation and velocity for two powder metal alloys IN-100 and Ti-6246 compared to 
a common engine disc material of nonpowder fabrication, i.e., Waspalloy. 


Table 2. 


Sample 

Velocity 
(mm/ pm) 

V L V S 

Grain 

Size 

(pm) 

Atten. 
(dB/ cm) 

“L “S 

Freq. 

(MHz) 

Freq. 

Dep. 

a s 

Waspaloy 

6.04 

3.23 

100 

4 

11 

10 

f 2 ' 2 

f 2 * 5 

IN-100 

6.52 

3.29 

2 

4 ± 1 

20 

50 

f 2 

£ 2 • 2 

Ti-6246 

6.11 

3.16 

2-4 

10 

- 

100 


^ 3 . 5 
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ANALYSIS 


This section presents the analytical developments leading up to the calcu- 
lations for a comparison between the experimental observations and the predicted 
frequency dependencies of back-scattering. 

The model for the microstructural backscatter (noise) spectrum we will use is 
taken from the work of Fertig and Richardson [ref. 9]. This model accounts for the 
statistical nature of the scattering process as well as the diffraction associated 
with the measurement system. Consider figure 11 below. We assume that the time 
window, T, associated with the experiment is small enough so that the variation in 
the sound beam along the z-axis is negligible within the sampled part. If we con- 
sider a scatterer at position r. = (r , , z,), with _r j ■ (x. ,y.), then the contribu- 
tion to the Fourier transform of! n(t), the measured J time dependent wave form, from 
this scatterer will be 


p(to,x^;z 1 ) exp(iodTj) Aj(io) 

where t- = (z. - zi)/(2c) is the relative time delay of the jth scatterer, A.(cj) is 
its scattering amplitude and p(o),_r.) (suppressing the dependence of zi) is tne 
system response function. 
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Figure 11 Diagram of pulse echo experiment measuring the 

backscatter from a region a distance from the 
transducer. 
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The spectral density of n(t) is found by taking the expectation of n(a))n*(oj) 
over the distribution of scatterer sizes and positions. If we assume that the scat- 
terers are uniformly distributed in real space, that multiple scattering is negli- 
gible, and we restrict the frequency mesh to = 2nk/T, we obtain 

oo 

E(n(w)n*(u)' )) = 6^, (JJ|p(u),_r) | 2 d 2 jr) p v e|a(w)| 2 , (6) 

— OO 

where p v is the expected number of scatterers per unit volume, 6 , is the Kronecker 

delta, and e|a(u))| 2 is the expected squared scattering amplitude^or a randomly 
chosen scatterer. For piston-source transducers, the integral of |p(w,_r)| 2 is not 
analytically tractable. For Gaussian transducers, however, it is quite easy. In 
either case, we identify the "volume" , vo1((jl>), of the beam, as a frequency dependent 
entity which may be defined by 

oo 

vol(to) = (cT/2) //|p(o),£) | 2 d 2 £ / 1 p(to, o) | 2 . (7) 

— OO 

Fertig and Richardson [ref. 9] give approximate expressions for vol(u)) for which 
oblique incidence and curved surfaces are allowed. The term p(u),o) is the on-axis 
diffraction correction and may be determined from the work of Thompson and Grey 
[ref. 10, 11] for piston sources. These latter authors give explicit forms for 
p (oo , o ) which allow for the inclusion of a reference calibration wave, R(u)), which 
enables one to relate a pressure profile on the face of the transducer to a voltage 
reading. In the case in which the scatterers are spherical and for which the 
Rayleigh limit holds, one obtains 

E | n(w) | 2 / |R(w) | 2 = 


2 C 2 ( Zl 

v 2 D(z i 


,z 2 ,u) 
»Z 2 


~k ~ 2 vol(w) 0 ) 4 rip E(r 6 ) 
R 1 R 2 


( 8 ) 


Here C and D are diffraction correction functions for the experiment and calibra- 
tion, Ri and R 2 are corner reflection coefficients associated with the calibration, 
E(r 6 ) is the expected value of the sixth power of the pore radius, and t] is a 
frequency and size independent parameter related to the scattering of spheres in the 
host material. In the far field limit, the expression in brackets above is 
independent of frequency. 


COMPARISON BETWEEN THEORY AND EXPERIMENT 


The experimental approach centers around a broadband transducer which is used 
to launch a short pulse of shear wave energy into the material in many different 
directions as shown schematically for the spin-pit test sample in figure 12. During 
its passage through the material, the microstructural inhomogeneities reflect energy 
back toward the transducer which is now used as a receiver to detect these small 
echoes which form a train of waveforms which have amplitudes and phase relationships 
depending on where and how strong the sources are. Several time traces recorded for 
a variety of regions for both samples A and C are shown in figure 13. It is clear 
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Figure 12 Schematic diagram of spin pit test specimen and ultrasonic 
shear wave beam. 

that these signals are more or less random and form a background. In a typical run, 
the backscattered background signals are windowed at 10 different sections of the 
time trace corresponding to 10 different regions of the sample. Each waveform is 
time averaged over 512 transmission pulses to decrease electronic noise in the re- 
ceiver equipment. The 10 waveforms are deconvolved with a reference signal to 
remove the transducer characteristics and averaged to produce one statistically 
representative waveform for the sample material from which now the power frequency 
spectrum may be calculated. Figure 14 shows the backscatter power spectra for 
samples A and C before deconvolution whereas figure 15 shows the same data after 
deconvolution and therefore presents the results in final form. 

As indicated by the dashed lines in figure 15 the predicted backscatter spectra 
have slopes which agree well with those obsrved, i.e., a 3.25 power law for the fre- 
quency dependence of the backscattered intensity. 

Another, more stringent, test of the data and model is to see whether a key 
microstructural feature may be predicted and compared with micrographical data. 

This is attempted in Table 3. Here the values of the estimated mean power densities 
for the two samples and the ultrasonic backscatter intensities at 10 MHz are in- 
serted into the model to calculate values of <r^>^ 6 which are then compared to 
corresponding values estimated from the micrographical examination. The comparisons 
are shown in the last two columns of the table. Considering that the micrographical 
examination was limited to the surfaces of the samples, the agreement is reasonable. 
The calculation demonstrated also that small errors in the count of the large pores 
introduces large errors in the final estimates. 
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Figure 13 Ultrasonic backscattering noise in IN-100. 
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Figure 14 Ultrasonic backscatter power spectra for two samples 
of IN-100 before deconvolution. 
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Pore Density 
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nr 
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(pm) 
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(-50 dB) 
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89 
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(-73 dB) 

5.01 x 
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FORMAL 

APPROACH TO 

INVERSE 

PROBLEM 



Whereas we have not attempted to solve the full inverse problem in this paper, 
it is appropriate to consider how one might attack this problem. For this purpose, 
it is convenient to formulate the problem in terras of a new set of random variables, 
called v .. , which represent the number of spheres whose centers lie in a volume 
element oVj located at r. in the material and who have radii in the range(a^ ,a ]• 
If we let ^(a^) represent the "counting” process equaling the actual number of 
spheres perunit volume with centers in 6Vj and radii less than or equal to a^, we 
see 
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jk 
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Again, if we take the case of negligible multiple scattering, the Fourier transform 
of the time windowed measured waveform takes the form 


y(w) = E E p(w,r j) A((o,a k ) v^ k 
j k 


+ e 


where e is additive noise (e.g., electronic), and ACtJ,^) is the scattering ampli- 
tude of a sphere with nominal radius a k e (a k _.,a ). The term p(to,r.) is the system 
response function for position r., the center 11 of 6Vj. The above formulation 
assumes that the 6Vj and 6a k = a k - a k _^ are a ll suitably small. 

Contrary to the "usual" inversion problem wherein information on an isolated 
scatterer of interest is contained in the expected value of y(w), for this case, the 
information on the statistics of the v^ k is contained in the variance of y(w). 

We make the following assumptions for the observation time T, 


1 . 

2 . 

3. 


p(to,r) = p(w,r) e i(JT , 


i> = 2z/c. 


~ stationary Poisson process 
e(o)) ~ zero mean stationary process, independent of v^ k . 


Then as 6V , o, Aa = max|a, - a, | o, we have on the special frequency mesh 
w k = 2itk/4, 


E(y(w k >) = o 

E(y(oj k )y*(u) k ,)) = 6 kk , [rjr j d 2 r |p(to k ,r)| 2 

. / da p(a) | A(to,a) | 2 + C £ (k)T] , 

where 6 kk « is the Kronecker delta, p(a) is the number density of spheres of radius 
a, and C (k)T is the variance of the k' th frequency component of the additive 
noise. e 

By inversion, we mean based on a set of samples of y(u)), estimate p(a). 
Specifically, we are interested here in estimating the underlying statistics of the 
pore distribution, not the actual size and location of a specific set of pores in a 
particular region of the specimen. This latter question is, of course, of interest 
in many applications. Here, however, we are interested in the former. As is well 
known, this is typically an ill-posed problem. Large changes in p(a) can result in 
small changes in y(w). For that reason, it will probably be necessary to 
parameterize p. To that end, define p-j as the number density of spheres 
and ¥ 3 (a;£) as the fraction of spheres with radii greater than a. Thus 
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p F (a;_0) “ / p(a)da . 
a 

The subscript 3 is meant to reinforce the fact that we are estimating volume 
densities of particles and not area densities. The latter often forms the raw data 
from sectioning experiments. The inversion problem is to estimate p 3 and the 
parameter 0_ (possibly vector-valued). One approximation to use for estimation of p3 
and 0_ is to assume that y((j) is a Gaussian process. This will be the case if the 
spheres are "dense" relative to the system response function so that no scatterer is 
"isolated" in time. In the Gaussian case, the posterior probability of p3 
and 0. given y is 


In P(p 3 ,£|y) a K - y*(w k ) B kk y(u> k ) 

k 


+ In P(p 3 ,i) +j Z In B kk , 

1c 

where K is a normalization constant, P(p3,0_) is the prior distribution of p3 and 
0., and is the inverse of the diagonal of the variance of yCw^)* It is given by 

B kk = T W + / P 3 f 3 (a;e)|A(u) k ,a)| 2 da , 


where 


M(w) = cT/2 / |p(u,r)| 2 d 2 j: 


and f 3 (a; 0 ) = -dF 3 (a;6^)/da is the probability density function of the sphere sizes. 

Two classical sets of estimators present themselves. The first set is |p3,0_} 
which maximizes the posterior density P(p 3 ,^|y). These estimators define the mode 
of this distribution and represent the most probable values. The second set esti- 
mators is {p 3 , 9j defined by the means of marginals of P(p3,0jy). Thus 

A 

A 

p 3 ■ / p 3 P( p 3 »ily) dp 3 d£ 


1-/0 p(p 3 ,ejy)dp 3 d0 


These estimators minimize the mean square error. 

In either case, the posterior variance of the estimators will be reduced if 
multiple waveforms y r (w), r = 1, ..., R are available. Having a small number of 
waveforms places a greater reliance on the prior density P(p 3 ,jO to keep the problem 
well-posed. 
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CONCLUSION 


The microstructure of a nickel-base powder metal alloy IN-100 has been studied 
with the aid of micrographical and ultrasonic techniques. The material showed 
grains of small diameter with a narrow distribution in the range of 2-4 pm and 
sparsely distributed micropores with diameter in the 5-20 pm range. The micropores 
dominated the ultrasonic attenuation and backscattering characteristics. The 
ultrasonic backscattering was selected to carry out quantitative measurements and 
calculations for comparison of theory and experiment. The frequency dependence from 
8 to 18 MHz of the backscattered signal was in good agreement with calculations 
which included the contributions from on-axis diffraction, back-surface reflection, 
the frequency dependent beam, the Rayleigh scattering and the distribution of scat- 
terers. In addition to the forward problem, the inverse problem was addressed in 
two ways (1) by estimating the mean pore size <r 6 > 1/6 from an approximate count of 
the pore density. This was found to give reasonably good agreement with the pore 
size estimated from micrographical information. (2) A formal approach for a 
stochastic inversion was developed. 

The results point the way toward the nondestructive detection and character- 
ization of anomalous distributions of micropores where conventional ultrasonic imag- 
ing is difficult. Such information may be useful in the early detection of poten- 
tial failure sites where stress induced void coalescence leads to crack initiation. 
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